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ABSTRACT 

The results of three-dimensional numerical simulations of mass transfer in semide- 
tached binaries after mass transfer termination is presented and the structure of 
residual accretion disk was investigated. It was shown that after the mass transfer 
termination the quasi-elliptic accretion disk becomes circular. It was also shown that 
the cessation of mass transfer results in changing of the structure of the accretion 
disk: the second arm of spiral shock is appeared as well as a new dense formation 
(blob), the latter moving through the disk with variable velocity. The blob doesn't 
smear out under the action of dissipation but is sustained by interaction with arms 
of spiral shock for practically all the lifetime of the disk. 

We also analyze the dependence of disk's lifetime on the value of viscosity. For 
the value of parameter a that is typical for observing accretion disk (a ~ 0.01) the 
lifetime of residual disk is found to be equal to 50 orbital periods. 



traditional 'hot spot' model. Successful appli- 
cation of the 'hot line' model for analysis of 
light curves in cataclysmic variables_(Bisikalo 
et al. 1998ai; Khruzina et al. 200113) allows 
us to conclude that the model appears to be 
adequate for considered systems. 

On the other hand, some observations show 
that for a large number of semidetached bi- 
naries there exist both epochs of constant 
rate of mass transfer and epochs in which 
the rate of mass transfer changes j-consider- 
ably Isee, e.g., Bath et al. 1924^; Wood 
19770; Bath & J^ringle 19820; £ath & 
van Paradijs 19830; Gilliland 19850; Rit- 
ter 1988[Ei_, Murray, Warner & Wickramas- 
ingheJiOOOlla; Schreiber, Gansicke & Hessman 
2000[cj). Evolutionary scenarios also predict 
the possibility of the contraction of the mass- 



1 Introduction 



Semidetached binaries are known to be inter- 
acting binary stars in which one component 
fills its Roche lobe and looses mass through 
the vicinity of inner Lagrangian point Li. Ear- 
lier we developed the 3D gasdynamical model 
for simulation of the mass transfer in semicLe- 
tachecL binaries.^ (see Bisikalo et al. 19970, 
1998b0, 1998ci). When applying to binaries 
with constant mass transfer rate this model 
allows to determine the main features of flow 
structure (see Bisikalo £t al. 19970, 1998b0, 
1998ci, 1999i, 200013). In particular, the 
stream interaction is shown to be shock-free 
for steady-state regime of the flow and the 
model of 'hot line' was introduced for bina- 
ries with constant rate of mass transfer instead 
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losing star. In this case the star doesn't fill the 
Roche lobe anymore and mass transfer will be 
terminated. 

The purpose of this paper is to investigate 
the flow structure in a semidetached binary 
after the mass transfer termination. On the 
first stage we have conducted 3D gasdynamical 
simulation of mass transfer for the case of con- 
stant rate of mass transfer up to the steady- 
state solution. Then we adopt that matter 
outflow is ceased and consider the structure 
of residual disk. It is clear that the evolution 
of accretion disk is controlled by physical pro- 
cesses responsible for the redistribution of the 
angular momentum in the disk. To investi- 
gate the influence of viscosity we conduct 3 
runs for various values of viscosity, these ones 
corresponding to following values of parame- 
ter a (in terms of the a-disk): a ~ 0.08 0.1, 
a - 0.04 ^ 0.06, and a - 0.01 ^ 0.02. 



2 Statement of the 
problem 

Let us consider a semidetached binary system 
and adopt that accretor has the mass Mi, the 
mass-donating star has the mass M2 , the sep- 
aration of the binary system is A, and an- 
gular velocity of orbital rotation is O. The 
flow of matter in this system can be described 
by Euler equations with equation of state for 
ideal gas P = {'j — l)pe, where P ~ pres- 
sure, p - density, e - specific internal energy, 
7 - adiabatic index. To mimic the radia- 
tive loss of energy we adopt the value of 7 
close to 1: 7 = 1.01, which corresponds to 
the near-isothermic case (Sawada, Matsuda & 
Hachisu 1986JE3 ; Molteni, Belvedere & Lan- 
zafame 1991 0; Bisikalo et al. 1997i). 

To obtain the numerical solution of the 
system of equations we used the Roe-Osher 
TVD scheme OjLa high approximation or- 
ffi,oe 19860; Chakravarthy & Osher 
with Einfeldt modification (Einfeldt 
The original system of equations 
were written in a dimensionless form. To 
do this, the spatial variables were normalized 
to the distance between the components A, 
the time variables were normalized to the re- 
ciprocal angular velocity of the system 
and the density was normalized to its value 
in the inner Lagrangian point Li. The gas 
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flow was simulated over a parallelepipedon 
[ V2A . . . X [- 1/2 A . . . 1/2 A] X [0 . . . V4A] 
(calculations were conducted only in the top 
half-space). The sphere with a radius of 
Vioo^ representing the accretor was cut out 
of the calculation domain. 

The boundary and the initial conditions 
were determined as follows: (i) we adopted 
free-outflow conditions at the accretor and at 
the outer boundary of the calculation domain; 
(ii) on the flrst stage in gridpoint correspond- 
ing to Li we injected the matter with param- 
eters p = p{Li), Vx = c(Li), Vy ^ ^ 0, 
where c{Li) is a gas speed of sound in Li 
point; (iii) for this stage we used rarefied back- 
ground gas with the following parameters po — 
10-5. p(Li), Fo = 10-V(ii)c'(J^i)/7, = 
as the initial conditions; (iv) on the second 
stage when steady-state regime is reached at 
the moment of time t = to we decreased the 
density of the injected matter to the value 
P = Po- 

The analysis of considered problem shows 
that the gas dynamical solution for the semide- 
tached binary is defined by three dimensicm- 
less parameters (Bisikalo et. al. 1998bO, 
19998; Lubow & Shu 1975y): the mass ra- 
tio q — M2/M1, the Lubow-Shu parameter 
e = c{Li)/An (Lubow & Shu igTSM), and 
the adiabatic index 7. The value of the adia- 
batic index was discussed above and we used 
the value 7 — 1.01. Analysis of J3ur presti- 
ous results (Bisikalo et al. 1998bi, 19990) 
shows that the main characteristic features of 
3D gas dynamical fiow structure are qualita- 
tively the same in wide range of parameters q 
and e. Therefore for the model simulation we 
chose them as follows: q = 1, e = Viq. 

To evaluate the infiuence of the viscosity 
on the solution, several runs with different 
spatial resolution were conducted. The Euler 
equations do not include physical viscosity so 
we varied the numerical viscosity by virtue of 
changing of computational grid. Three grids 
were chosen for simulation: 31 x 31 x 17, 
61 X 61 X 17, 91 X 91 X 25 (hereinafter runs 
"A", "B", and "C", correspondingly). In 
terms of the a-disk the numerical viscosity 
for runs "A", "B", "C" approximately corre- 
sponds a ~ 0.08 ^ 0.1, a - 0.04 ~ 0.06, and 
a - 0.01-^0.02. 
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Fi sure 1 '. Time evolution of the total mass of the 
calculation domain (top part of the figure) and the 
mass of the disk (bottom part of the figure) for "A", 
"B", and "C" runs. The mass of uniform matter with 
background density p = po for the whole computation 
domain and for the disk are shown by dashed lines. 

3 Flow structure after the 
mass transfer 
termination 

Prior to simulation of the flow structure with 
'turned-off' mass transfer the near-steady- 
state solutions for the case of constant non- 
zero rate of mass tiansfer were obtained 
(Bisikalo et al. 2000llj) and used as initial 
conditions. At the moment of time t = to the 
rate of mass transfer was decreased in five or- 
der of magnitude, which is correspond to the 
cessation of mass transfer. 

Calculations of runs "A" and "B" were 
lasted until the density of gas becomes less 
than the background density po = 10~^ p{Li). 
This corresponds to the full vanishing of mat- 
ter due to accretion and outflow through outer 
boundary. The durations of these two runs 
correspond to 5Porb and 12Porb after the mo- 
ment of time tg. Run "C" has the best resolu- 
tion and minimal numerical viscosity. Nev- 
ertheless this run was conducted to lOPorb 
only, since this run is very computer-time- 
consuming while the time when disk vanishes 
can be estimated as ^ 50Porb- We extrapo- 
lated the results of calculation of run "C" for 

t > lOPorb. 

Before investigating of the structure of 



gaseous flows in binary system let us estimate 
the lifetime of the residual disk after the termi- 
nation of mass transfer for different values of 
viscosity. The lifetime of the residual disk can 
be evaluated from the Fig. 1 where the total 
mass of the calculation domain (top part of the 
figure) and the mass of the disk (bottom part 
of the figure) versus time is presented. The 
lower limit of the density of gas is the back- 
ground density pa, therefore the mass of the 
uniform matter with p = pQ calculated for the 
whole computation domain and for the disk 
are asymptotic lines for graphs M{t). These 
lines are shown in the Fig. 1 by dashed lines. 

The analysis of the data presented in the 
Fig. 1 shows (and there is no surprise now) 
that the lifetime of the residual disk is in- 
creased when the numerical viscosity is de- 
creased: for run "A" {a ^ 0.1) Tdisk — 5Porb', 
for run "B" {a ^ 0.05) Tdisk = 12Porb] and for 
run "C" (a ~ 0.05) the extrapolation of the re- 
sult of simulation shows that the lifetime of the 
residual disk exceeds 50Porb- It is pertinent to 
note that the value a ~ 0.01 is typical for ob- 
servable accreticm,disks (see, e.g., Lynden-Bell 
& Pringk^l9740; Meyer-Hofmeistei_& Rit- 
ter 19330; Armitage & Livio 1996§; Tout 
1996[l3). Consequently, we can expect the 
lifetime of residual disk after the termination 
of mass is of order BOPorb (near a week) for 
typical dwarf novae. 

Now let us consider the evolution of the 
residual accretion disk in time. Since the re- 
sults of runs "A", "B", "C" are quahtatively 
similar, we will focus on the run "B" and ana- 
lyze it in detail. The choice of run " B" is called 
forth by the following reasons: firstly, the spa- 
tial resolution of this run is well enough to 
catch all the details of the disk structure, and, 
secondly, the lifetime of residual accretion disk 
is sufhciently small {12Porb) so we can study 
its fate up to full vanishing. The Figure 2 
presents isolines and bird-eye views of density 
in the equatorial plane for different moments 
of time, the upper panels correspond to the 
moment of time t — tQ + 0.08Porb, middle pan- 
els - to t — to + O.ldPorb, and lower panels 
- to t = to + 0.9Porb- The analysis of these 
figures shows that for initial moments of time 
{t < to + O.OSPorb) the structure of accretion 
disk doesn't differ much from the structure ob- 
tained for the steady-state solution with con- 
stant mass transfer rate (see Bisikalo et al. 
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Figure 2: Isolines (left panels) and bird-eye vies (right panels) of density in the equatorial plane. Three rows 
of panels correspond to the following moments of time: t = to + Q.OSPorb, t = to + 0.16Pori)i t = to + O.QPort, 
where to is the moment of the mass transfer termination. 
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Figure 3: Time variation of the mean density of mat- 
ter passing through a semi-plane XZ {Y = 0, X > A) 
shcing the disk. A zoom for the mean density variation 
for times from io to -I- lAP^rh is shown in the top of 
the figure. Two sets of moments of time (8 moments 
per set, each set covers one period of the mean density 
variation) are also shown. 

1997!, 1998bS, 1998ci, 1999!, 20000). 
As always the stream of matter from Li domi- 
nates. It is clearly seen that the uniform mor- 
phology of the system stream-disk results in 
quasi-elliptical shape of the disk and the ab- 
sence of 'hot spot' in zone of stream-disk in- 
teraction. At the same time the interaction 
of the gas of circumbinary envelope with the 
stream results in the formation of an extended 
shock wave located along the stream edge ('hot 
line'). The Figure 2 also manifests the forma- 
tion of tidally induced spiral shock located in I 
and II quadrants of coordinate plane. Appear- 
ance of tidally induced two-armed spiral shocik 
was diacoveredjJjL Sawada et al. (1986alc3, 
1986bEl, 19870). Here we see only the one- 
armed spiral shock. The flow structure in the 
place where the second arm should be is de- 
fined by the stream from Li which appears to 
suppress the formation of another arm of the 
tidally induced spiral shock. 

From the panels of middle row of the Fig. 2 
it is seen that after the moment of time t ~ 
tf) + Q.lQPorb the stream from Li has no in- 
fluence on the structure of residual accretion 
disk anymore and the disk is circularized. It 
is also seen that the shock wave located along 
the stream edge ('hot line') disappears and the 
second arm of the tidally induced spiral shock 
located in III and IV quadrants of coordinate 



plane is formed. It is illustrated by lower row 
of the Fig. 2 where results of simulations for 
the moment of time t = -\- 0.16Por& are pre- 
sented. Note that the intensity of shocks as 
well as the mass of accretion disk decreased in 
time but the qualitative flow structure retains. 
The flow structure is changed when the tidal 
interaction doesn't result in the formation of 
spiral shocks anymore (it corresponds to the 
moments of time + lOPorb + H-Porb for 
run "B" ). Therefore, we can conclude that the 
spiral shocks exists during of all lifetime of the 
disk (remind that it is of order VlPorb for run 
"B"). 

Our results also show the formation of dense 
blob in the residual disk. The blob is formed 
after the stream vanishes, i.e. after time 
~ to -|- 0.15Por-6 and rotates with variable ve- 
locity. The Figure 3 presents time variation of 
the mean density of matter passing through a 
semi-plane XZ (Y = Q, X > A) slicing the 
disk. Along with a general density fall we can 
see here periodic oscillations called forth by 
passing the blob. The blob doesn't smear out 
under the action of dissipation and the period 
of its rotation (^ O.lSPorb) remains the same 
practically up to the moment when the disk 
vanishes. Figure 3 also contains a zoom for 
the mean density variation for times from tf) to 
tQ+lAPorb- Here we marked 2 sets of moments 
of time (8 moments per set), each set covering 
one period of the mean density variation. The 
flrst set corresponds to the moments of time 
to + O.llPorb ^ ^0 + 0.26Porb, and the second 
- to to + 0.91Porb ^ to + IMPorb- Figurcs 4a 
and 4b show the density distributions for the 
moments of time of these two sets. Analysis of 
the data of Fig. 4a, b, shows that the blob re- 
tains by its interaction with the arms if spiral 
shocks. 

Let us consider this mechanism in de- 
tail. Initially appearing as the residue of the 
stream, the blob tends to distribute uniformly 
along the disk under the action of dissipa- 
tion but retards after passing trough spiral 
shock so the compactness of the blob retains. 
Later the increasing density contrast between 
the blob and the disk and growths of pres- 
sure gradient will force the standstilled blob to 
move. When reaching the second arm of spiral 
shock the process of formation/retainment of 
the blob is repeated. This mechanism is con- 
flrmed by data of Fig. 5. Figure 5a presents 
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Fi gure 4b! Density distribution in the equatorial plane for the second set of moments of time from Fig. 3. 
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the; velocity of the blob (the velocity of the 
point with maximal density) versus time (solid 
line). The keplerian velocity of this point 
Vk = \/ GMi/r is also shown by dashed line. 
The velocity variations are shown for time in- 
terval to+0.91Por6-^io + l-ll-Por6 including the 
time interval for Fig. 4b and covering the full 
period of blob's revolution. Figure 5b shows 
variation of maximal density in the equato- 
rial plane for the same time interval. It is 
seen from Fig. 5 that matter is retarded after 
passing the shock and becomes more dense, 
the maximal density pmax being increased 1.5 
times. When accumulating a sufficient amount 
of matter, the blob comes off the shock and its 
velocity increases while maximal density de- 
creases. After passing the second arm of spi- 
ral shock the blob is retarded again. Note that 
the contrast of density between the blob and 
the disk remains the same up to the moment 
when the blob disappears. In this moment of 
time the disk is contracted and spiral shock 
disappears as well. Later the density contrast 
blob/disk decreases under the action of dissi- 
pation but the blob exists up to the full van- 
ishing of the disk. This is illustrated in Fig. 6 
where time dependency of the density contrast 
blob/disk is shown. The large value of den- 
sity contrast as well as the variable velocity of 
the blob revolution (but with constant period 
O.lSPord) make this feature very interesting 
for observations. 



4 Conclusions 

We have presented the results of 3D numer- 
ical simulations of mass transfer in semide- 
tached binaries after the mass transfer ter- 
mination. Prior to simulation of flow struc- 
ture with 'turned-off ' mass transfer the near- 
steady-state solutions for the case of constant 
non-zero rate of mass transfer was obtained 
and used as the initial conditions. At the 
moment of time 1 = 1^ the rate of mass 
transfer was decreased in five order of mag- 
nitude, which corresponds to the cessation of 
mass transfer. To investigate the influence 
of the viscosity we conduct 3 runs for vari- 
ous values of viscosity, corresponding to fol- 
lowing values of parameter a (in terms of a- 
disk): a ~ 0.08 0.1, a ~ 0.04 ^ 0.06, and 
a ~ 0.01 ^0.02. 

The investigation of structure of the resid- 
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Figure 5: Upper panel: A solid line - variation of 

the velocity of the blob (the velocity of the point with 
maximal density). A dashed line — keplerian velocity 
of this point Vk = \J GM\/r. The velocity variations 
are shown for time interval to+O.^lParh'^to + ^-^^Porh 
including the time interval for Fig. 4b and covering the 
full period of blob's revolution. 

Lower panel: Variation of the maximal density of 
the blob. 
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Fi gure 6; The density contrast between the disk and 
the blob versus time. 



ual disk reveals an essential increasing of its 
lifetime when the viscosity deceases. Even for 
a 0.05 the lifetime of residual disk exceeds 
12Porb, and for a 0.01 (that is typical for ob- 
servable accretion disks) the lifetime of resid- 
ual disk is as much as SOPorb- 

Our simulations show that near the moment 
of time t = to + 0.2Porb the flow structure 
is changed significantly. The stream from Li 
doesn't dominate anymore, and the shape of 
accretion disk changes from quasi-elliptical to 
circular. The second arm of tidally induced 
spiral shock is formed while earlier (before 
the termination of mass transfer) it was sup- 
pressed by the stream from Li. The mass 
of the disk is gradually decreased but spiral 
shocks exist practically up to the full vanish- 
ing of the disk. 

Our simulations also show that a dense blob 
is formed in the residual disk, the velocity of 
the blob motion through the disk being vari- 
able. The blob doesn't smear out under the 
action of dissipation but is sustained by in- 
teraction with arms of spiral shock up to the 
moment of time when spiral shock disappears. 
The density contrast between the blob and the 
disk is rather large ~ 1.6 and begins to de- 
crease only when spiral shock disappears (it 
occurs at time to + lOPort) but even for this 
time and up to the full vanishing of the resid- 
ual disk the density contrast is of order 1.2. 
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